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' Matter waves originating from a localized region in space appear commonly in 

^ physics. Examples are photo-electrons, ballistic electrons in nanotechnology de- 

, vices (scanning-tunneling microscopy, quantum Hall effect), or atoms released from 

■ a coherent source (atom laser). We introduce the energy-dependent Green function 

' as a suitable tool to calculate the arising currents. For some systems experimental 

^1^' data is available and in excellent agreement with the presented results. 
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1 Introduction 



The propagation of matter waves in external fields shows a wealth of interest- 
ing quantum-mechanical phenomena. Here, we only discuss uniform, static 
gravitational, electric, and magnetic fields. The corresponding Hamiltonians 
are quadratic polynomials in the momenta and position operators. Therefore 
all classical equations of motion are readily available, as well as the time- 
0^ ' evolution operator in position spac^ 

o 

m. X(r,t|r',io) = (r|exp(-iH(i-to)A)|r')- (1) 

o . 

However, the knowledge of the quantum-mechanical time-evolution is not suf- 
^ • ficient to describe experiments in which the energy E of the travelling parti- 

cle is fixed, rather than the travel-time t — Iq. The energy-dependent coun- 
■ terpart of the time-evolution operator is the energy-dependent Green func- 

I tion, as we will see below. Surprisingly, analytic results for three-dimensional 

^ energy-dependent Green functions are scarce compared to the available time- 

dependent kernels: Commonly applied methods like the Feynman-path inte- 
gral cannot be used. 

If the process is to be treated time-independent or stationary, we have a 
}_j ■ system with a reservoir (or quantum source) that feeds particles into the exter- 

I nal fields at a constant rate. This concept might be surprising in a quantum- 

mechanical context, since the Schrodinger equation implies a conservation of 
the probability current. However we will see that in the quantum-source ap- 
proach the currents are conserved outside the emitting source region. This is 
just one advantage compared to another possible approach: On can also solve 
the Schrodinger equation for a bound particle in the presence of an exter- 
nal field that enables the particle to tunnel from the originally bound statJ^. 
Since the overall time dependence of a decaying state is given by exjp{—iEt/fij, 
a negative imaginary component of E will lead to an exponentially decreasing 
probability current. While this feature is not unexpected in dealing with a 



1 



decay-process, the resulting wave-function and probability current cannot be 
properly normalized. In the following we consider stationary currents which 
implies that we can use a real-valued energy-parameter. In connection with a 
quantum source we can define a meaningful probability current. Historically, 
the quantum-source approach was advocated by SchwingeJ^. 

2 Quantum Sources and Currents 

Here we merely state some consequences of the intr odu ction of an inhomoge- 
neous source term a(r) to the Schrodinger equatiorPEl 

[S-H]V'sc(r) = (T(r). (2) 

The simplest example of such a source term is the (5-distribution, which de- 
scribes matter waves originating from a point-source. A point source allows 
us to obtain a solution in the form of a Green- function G{r, r'; E) 

[E-H]G{r,r';E)=d{r-r'). (3) 

In order to define the Green-function uniquely, we have to specify boundary 
conditions. In the following we impose the boundary-condition that we de- 
scribe the emission of outgoing waves. The energy-dependent Green function 
is linked to the time-evolution operator by a Laplace transform: 

/.oc 

G(r,r';£;) = lim / dt {r\U{t,to)\r') e''^'/'''-'^'^'' . (4) 

We obtain the scattering wave function assigned to the quantum source a{r) 
by a convolution integral: 

Ac{r) = J dV'G(r,r';£;)a(r'). (5) 

A quantity of interest is the associated current, its spatial distribution and 
energy dependence. The current distribution is actually experimentally ac- 
cessible, and we will discuss some examples. We define the current density 
distribution j(r) 

j(r) = A3[^^^(r)*VVsc(r)] - ^^|^sc(r)|^ (6) 
TO m 

where A(r) denotes the vector potential. Integration of j(r) over a surface 
enclosing a{r) yields the total current J{E) emitted by the source: 



JiE) = -|3 



J dV J d^r'cr{r)*G{r,r';E)cr{r') 



(7) 



Connected to the appearance of a source for particles we have to modify the 
equation of continuity. Instead of V • j(r) = 0, valid for a stationary system 
in the absence of sources, we now find: 

V.j(r) = -^SKrrVsc(r)]. (8) 
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Note, that outside the source region the current is a conserved quantity. 



3 The twin-slit experiment 

3.1 Material slits 

As a first example we address the propagation of particles from a point source 
in a field- free environment. We want to calculate the current distribution for 
a particle which is emitted from the source located at ta and arrives on the 
detector screen at r^. In between, a wall with two holes at , r_f/2 is placed. 
To simplify our considerations, we will assume a point source of particles. 
Since the double-slit configuration contains two possible particle paths, one 
from — > r^f^ — > and the other one from r//^ tb, we have to 

form a superposition of the two transition amplitudes 

if3iit(rB,T|r^,0) = / dii i^ft.ec(rB,r|rH,,ti) iffrcc(rHi,ti|rA,0) (9) 
Jo 

+ / dt2 K{^cc{rB,T\rH2,t2) K{rcc{rH2,t2\rA,0). 
Jo 

By integrating over ti and ^2 from to T we consider a time evolution that 
starts at i = at and ends at i = T at r^. In between a passage through 
the holes occurs at some intermediate time ti. Knowing the form of the 
propagator Kf^cc, we have to carry out the intermediate time integration over 
ti. However, in an experimental setup we would like to eliminate the need to 
measure exactly at time T. Instead one prefers to record a static pattern on 
the detector screen for particles emitted with a fixed energy E. Switching to 
the energy-dependent Green function we obtain 

i f°° 

Gsiit(rA, rs; ^ J dT e-'^^/''if,iit(rA, T\rB,0) 

= --/ dTe-'^^/'*^ / dt, Ki,,,{rB,T\rH^,t,) Ki,,,{rH„t,\rA,0) 
Jo j^^i Jo 

(10) 

Using the convolution theorem for the Laplace transfornIS! we can rewrite the 
Green function as 

2 

Gsiit{rA,rB;E) = i?i^ Gfroo(rB, r//,; i;) GhccirHi,rA;E). (11) 

1=1 

The last form clearly shows the superposition of two scattering waves. The 
corresponding current density distribution shows the familiar interference pat- 
tern, which was observed in a twin-slit experiment with electron!^. 
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Figure 1. Trajectories in the uniform force field F. The dashed parabola denotes the 
classical accessible region. In this region, two trajectories with the same energy connect 
the source S with a point on the detector plane. For electrons emitted in an electric field 
macroscopically visible interference fringes emerge. 



3.2 A virtual twin-slit provided by a homogeneous force field 

In 1981 it was pointed out that a linear force field can act as a double slit 

on 

without the need for an actual material slit . An extended analysis in 
the language of quantum propagators is 

availabltPm. 

The time-dependent 
propagator for a particle in a uniform force field F = Fe^ is given by 

,, , /TO \3/2 /imlr-r'P iFT , „ iF^tA 
K.Mr.T\r',0) = (^) exp (^^^ + + - , 

(12) 

where T = t — t' . As before, we want to calculate the probability amplitude 
for a particle traveling from the source to hit the detector at r^- Using 
again a monochromatic particle source, the stationary probability amplitude 
is derived similarly to H1U|I and reads 

Gfieid(r,r';^) = -i^ dT e'^^/''Xfieid(r, T|r', 0). (13) 

Although Gficid(r, r'; -E) is available in closed analytic for nP, we will retain 
the integral form for the following discussion. We can approximate the inte- 
gral with the method of stationary phases for a wide range of parameters. At 
a stationary point the derivative with respect to T of the exponent in H13(l van- 
ishes and we get a large contribution from the region close to this point to the 
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Figure 2. Plot of the current density distribution for parallel electric and magnetic fields. 
Parameters: Emission energy: E = 60.8 /leV, electric field: F = 116 eV/m, magnetic field: 
B = 0.001 T. There is rotational symmetry about the 2-axis. 



otherwise oscillatory integral. The resulting biquadratic expression in T has 
two roots for classically allowed motion E > 0, and thus we obtain two real- 
valued solutions for the time of flight Ti. Galile^l noticed this in his studies 
of the parabolic motion of free falling masses. The two classical trajectories 
produce the pronounced interference fringes depicted in Fig. ^ The spacing 
of the fringes depends on the energy of the emitted particles E, the strength 
of the uniform force field F and the distance between source and detector 
(r, z). Following the idea of Demkov et aLl^iSl Blondel constructed a device to 
realize the field d ouble-slit (also for electrons) , the so-called photodetachment- 
microscop j ^ "^ 1 The photo-detachment is a high-precision spectrometer for 
the determination of the energy parameter and is used to determine the elec- 
tron affinity of negative ions with extreme accuracy. 

4 Uncertainty principle in parallel magnetic and electric fields 

The addition of a homogeneous magnetic field to an electric field leads to 
additional classical trajectories that connect the source with a given point on 
the detector A parallel magnetic field will enforces a cyclotron motion of 
the emitted electrons and establishes a lateral confinement (see Fig. EJ. 

For low energies the classically allowed paths are confined closely to the 
path leading directly downwards from the quantum-source. In quantum me- 
chanics, the focussing of the electrons due to the electromagnetic fields cannot 
be perfect, since Heisenberg's uncertainty principle enforces a minimum lateral 
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Figure 3. Transition from a strongly confined condensate to a more extended source 
distribution. There is rotational symmetry about the vertical axis. The source width 
is denoted by a. Interference fringes are clearly seen for a < 0.4 fim. Parameter: 
Au = E/(2-Kh) = 2.5 kHz, F = m^b g, with g = 9.81 m/s^, and m^b = 87 u. 



distribution in position space along the escape path. Also negative energies 
E < lead to a tunneling current for which no classical trajectory exists. 
While the exact solution in terms of the energy-dependent Green function is 
still available, also a heuristic description in terms of the minimum uncertainty 
principle is possibl^^ and can be used to establish a tunneling time. 



5 Coherent atomic ensembles and the atom laser 

The quantum-source approach is not limited to point-like emitters. Ana- 
lytic solu tion s are available for a Gaussian source in an external homogeneous 
force fielcJlil. A Bose-Einstein condensate (BEG) is a possible realization of 
a macroscopic quantum source. The controlle d and coherent release of atoms 
from such a condensate is called atom-laser. Ir^^^we establish a completely 
analytical model for the description of the efficiency of such an atom laser 
and also for the density distribution in the resulting atom beam. Especially 
for strongly confined Bose-Einstein condensates (realized in so-called micro- 
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Figure 4. The transition from a weak electric field to a stronger electric field is shown. 
Parameter: Magnetic field; B = 0.5 T, parallel electric field; = 0, perpendicular electric 
field: Fy = (1,100,400) eV/m (dotted, solid, and dashed line respectively. 



traps) the appearance of strong interference phenomena in the atomic beam 
is predicted. Remarkably, in the far-field sector in a uniform force field the 
Gaussian source may be replaced by a virtual point source which is shifted 
upwards in the external field. Thus we can readily apply the analysis of the 
virtual twin-slit in the previous section. However, the shift in position affects 
the energy-parameter of the Green-function and leads to effective negative 
energies. Therefore, for larger condensates the interference fringes are sup- 
pressed (see Fig.OJ. 

6 Electromagnetic fields and the density of states 

Another use of quantum sources is the field of solid-state physics. There, 
considerable interest exists in expressions for the density of states (DOS) of a 
system. The DOS is defined by 

n(r; E) = {r\S{E - H)|r) = --3[G(r, r; E)]. (14) 

TT 

In our approach the DOS is directly proportional to the total current gener- 
ated by a point-source, which is an surface-integral over the spatial current- 
distribution around the source. For the spatial current, classical trajectories 
can be compared with the quantum solution and show the difficult transition 
between a semi-classical and purely quantum- mechanical regime. In Fig. ^ we 
demonstrate the influence of external fields on the DOS. 
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7 Outlook and conclusions 



Quantum sources and the energy dependent Green functions provide an excel- 
lent tool for the analysis of stationary propagation in external fields. It is pos- 
sible to extend the theory to include sources with an angular momentum (i.e. 
rotating BECs, or electrons coming from a specific atomic orbital^. In solid 
state physics, we can use the source-formalism to des cribe the electron prop- 
agation in a two-dimensional quantum-Hall systerr^l These examples surely 
present only a glimpse of the possible applications of the source-frameworkl^. 
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